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Abstract 



We study interfacial behavior of a lamellar (stripe) phase coexisting with 
a disordered phase. Systematic analytical expansions are obtained for the 
interfacial profile in the vicinity of a tricritical point. They are characterized 
by a wide interfacial region involving a large number of lamellae. Our 
analytical results apply to systems with one dimensional symmetry in true 
thermodynamical equilibrium and are of relevance to metastable interfaces 
between lamellar and disordered phases in two and three dimensions. In 
addition, good agreement is found with numerical minimization schemes of 
the full free energy functional having the same one dimensional symmetry. 
The interfacial energy for the lamellar to disordered transition is obtained 
in accord with mean field scaling laws of tricritical points. 
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I. INTRODUCTION 



Some equilibrium phases of matter exhibit spatial modulations in one or more of 
their local properties (local order parameter) and hence are called modulated phases ]T[] . 
Modulated phases can be found in a variety of physical and chemical systems such as 
diblock copolymers, magnetic garnet films, ferrofluids, adsorbates on surfaces of metals, 
superconductors (Type I) and semiconductors. Although these systems all have spatial 
modulations in their order parameter, the microscopic origin of these modulations is very 
different. 

Usually modulated phases manifest either a one dimensional packing of sheets (the so- 
called lamellar or stripe phase), two-dimensional packing of cylinders with an hexagonal 
symmetry or a cubic packing of spheres or other, more complex, bi-continuous super- 
structures with cubic symmetry. The origin of the preferred length scale of the modula- 
tions can be explained in many cases by a competition between two type of interactions: 
a short-range interaction which tends to make the system more homogeneous together 
with a long-range one, or a non-local one, which prefers proliferation of domain walls. For 
example, in the magnetic films, the dipole-dipole interaction has the latter effect. 

These two opposing interactions are responsible in the above systems for the appear- 
ance of spatially modulated phases having an optimal non-zero wave-vector. Depending 
on the range of the thermodynamical variables, instead of having the usual segregation of 
the system into two homogeneous phases (such as coexisting liquid and gas phases, or two 
coexisting polymer phases), the system self-organizes into alternating domains. Those 
domains are associated with spatially modulated order parameter (polarization, magne- 
tization, relative concentration), and can be regarded as unit cells of super-structures of 
lamellar, hexagonal or cubic symmetry. Depending on the system in mind, the domain 
length scale can be as small as one hundred angstroms or as large as centimeters [[j]. In 
contrast to periodic crystalline structures, the periodicity in modulated phases can be 
tuned by varying external fields such as pressure, temperature and electric or magnetic 
fields. 

In this paper we would like to present analytical results related to interfaces in mod- 
ulated phases. When the modulations are relatively weak an analytical expansion of the 
interfacial profile can be obtained systematically. In particular, we calculate the profile 
of the lamellar-to-disordered interface for a model free energy where spatial modulations 
are restricted along one dimension. Physical systems can create more complicated modu- 
lations in two- and three dimensions. In many cases, the lamellar-to-disordered interface 
addressed here can be shown to be only metastable and not in true thermodynamical 
equilibrium. Nevertheless, their systematic study can be useful to understand nucleation 
and growth in e.g., polymeric systems with slow dynamics. In addition, the rather sim- 
plified lamellar-disordered interface can be generalized to other situations. 

Our analytical results are valid in the vicinity of a tricritical point which occurs for 
asymmetric profiles in our strictly one- dimensional model. Close to the tricritical point it 
is possible to decouple the two characteristic length scales which appears in the model : 
the period associated with the modulated phase and the width of the interface. Thus, it 
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is possible to compute analytically the leading terms of the interfacial energies. 

The obtained results presented here are in good agreement with recent work by Netz 
et al |2| where numerical minimization of a similar free energy functional yielded detailed 
structure of interfaces between various two-dimensional modulated phases. 

II. THE MODEL AND ITS PHASE DIAGRAM 

A phenomenological free energy F similar to the one introduced in Ref. |2| is used 
hereafter in order to allow us to study one- dimensional profiles between lamellar and 
disordered phases. 

F = 2(V 2 0) 2 - 2(V0) 2 + \<? + \^ (1) 

This mean-field free energy is a modified Ginzburg-Landau (GL) expansion in the order 
parameter 4>(x), which varies between = 1 for A-rich domains to <fi = for B-rich 
domains. The last two terms in F are the usual terms appearing in the Landau expansion 
with t = (T — T c )/T c being the reduced temperature and the coefficient of the fourth order 
term is chosen to be 1/4 without loss of generality. For an order parameter without spatial 
modulations they define our disorder free energy, F D (within the Landau expansion): 

F D = + (2) 

The main deviation from the regular GL expansion can be seen in the first two terms of ([!]). 
The coefficient of the gradient square term is negative (favoring modulations), whereas the 
coefficient of the Laplacian square term is positive (preventing strong modulations) and 
playing the role of a homogenizing interaction. An illustrative example can be a diblock 
copolymer melt where each polymer is composed of two blocks, A and B, which would 
like to phase separate. However, since the A and B chains are chemically linked forming 
one polymer chain, the system can only undergo a micro-phase separation resulting in a 
spatial modulation in the local A/B relative concentration. 

In order to construct the phase diagram, we compare the mean-field free energy F D 
associated with the disordered phase (0 = const.) with the free energy F L of a phase 
described by a modulating <p{r). Our approach involves several simplifying assumptions: 
In the vicinity of a critical point (weak segregation limit), an analytical expansion of the 
free energy in powers of the small order parameter (Landau expansion) can be justified. 
Since a mean-field approximation is employed, critical fluctuations corrections are omitted 
H|J. Furthermore, only the most dominant modulation mode (the so-called single-mode 
approximation) is considered because the amplitude of the fluctuations of other g-modes 
grows more slowly in this region as was shown by Brazovskii Finally, only one- 

dimensional modulations of 0(r) = <p(x) along an arbitrary chosen x direction will be 
explored. This reduces the problem into an effective one-dimensional problem. It applies 



3 



to interfaces where the lamellae are themselves parallel to the lamellar - disorder interface. 
We further discuss the validity and limitations of our assumptions in the last section. 

4> = 4> + 2<p q cos(qx) (3) 

where O is the mean spatial order parameter (the average A/B relative concentration) 
and 4> q is the amplitude of the g-mode modulation. 

The most dominant mode is obtained by minimizing the free energy with respect to 
the q- vector. Within the single mode approximation this leads to q* = 1/ y/2 with our 
choice of simple coefficients in ([!]). This is the wave- vector characterizing the periodicity 
of the super-structure, 27r/g*. Inserting <j>(x) of the dominant mode from (|3]) into ([[]) we 
obtain the lamellar free energy 

^ = ^+^+(^-1+3^+^: (4) 

The disorder free energy F D (0) is simply given by setting <fr 2 = in the above expression 

lor /•; . 

Since this expression is symmetric with respect to <ft — > — <p only the O > region 
can be considered without loss of generality. We then apply the variational principle with 
respect to <fi q . This enables us to integrate over <p q — one of the two order parameters. 
Since the amplitude of the modulation minimizing F L is : 

the necessary condition for the existence of this modulated phase is <fi 2 > 0, hence 

^-#>0 ( 6 ) 
Substituting <p 2 from (^|) into (^) gives the free energy for the lamellar phase 

f l = + (i - \w - yi <7) 

In the region bounded by inequality @ and close to the critical point, the free energy 
is a convex function of <p . Therefore, the second-order phase transition line separating 
the order and the disordered phases is given by: 

4>l = ^ (8) 

showing an upward shift of the critical temperature t c = 1 from its usual value for the GL 
model, t c = 0. A tricritical point appears at a temperature t tcp = §, for which 5 2 F L /8(j) 2 l = 
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0. Below this temperature, the free energy is concave; the system therefore prefers to 
phase separate. Equation gives the value of the order parameter (concentration) at 
the tricritical point, 4> 2 T = 2 (tj cp ) = ^. 

In Figs. 1 and 2 the phase diagram is plotted as function of the reduced temperature 
- average concentration (t, O ), and as function of the chemical potential - reduced 
temperature (/x, t), respectively. The phase diagram is calculated numerically within the 
single mode approximation by comparing the lamellar free energy F L with the disorder 
one, F D . We note that the phase diagram in the vicinity of the tricritical point resembles 
that of the BEG (Blume- Emery- Griffiths) spin one model 0. The BEG model also has 
two order parameters and exhibits a similar behavior: a first order phase transition line 
meets a second-order one at a tricritical point. 

In the vicinity of the tricritical point (but inside the two phase region) we are still in 
the weak segregation limit. Hence, a systematic expansion can be performed as function 
of a new reduced temperature e, measuring the distance from the tricritical point: e = 
ttcp — t > 0. We then get : 

f l = 1<pI + \4-1(€ + 1-<pD 2 0) 

This equation can be expanded in terms of r\ D and r) L , such that (fi(e) = (p T + r) D (e) for 
the disordered (isotropic) phase and 0(e) = T + r] L (e) for the lamellar phase. In the 
two-phase coexistence region, the boundaries are given by the common tangent method : 

V D = Vl = V ( 10 ) 

together with 

G l = F l ~ = G D =F D - fj4> D (11) 
where F L and F D are expressed in term of r] D and r] L , respectively. 

F L = ~/ T + + + MtV* + ~VI - 50 T ^ - %f L (12) 

F D = ~/ T + ~/ T + {\~ e)<P T V D + (1 - e)i + <j> T r,l + \rf D (13) 

Relation (|T0D enables us to write r] D as an expansion in r\ h up to third order in e : 
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and relation ( 



IIP leads to r] L up to third order in e 



Therefore 



V D =4> T (2s + -/-^) (16) 
In terms of newly defined order parameters: <fi± = ((f) D ± <p L )/2, we get 

0- = ^ D - 0J = \(v D - V L ) = * r (f e + ^ - (17) 

n 1 £ ,2 2 6 3o 

C = <jt = -e e d 

r9 3 5 50 500 

Note that at the tricritical point, the slopes of the two phase transition lines given by 
Eq. (|8|) and Eq. (p~6|) are equal as has been noticed by the original mean- field work of Ref . 



III. THE LAMELLAR-DISORDER INTERFACIAL PROFILE 

In a simple condensation transition, the interface width between the two homogeneous 
phases, £ ~ t -1 / 2 , has the same scaling behavior as the typical concentration fluctuations 
of the homogeneous phase. In a modulated phase, a second characteristic length scale 
exists and has to be considered together with this correlation length. This is the wave- 
length of the super structure, q* . In order to compute analytically the profile of the order 
parameter for the disorder- lamellar interface close to the tricritical point, where the two 
length scales are expected to be uncorrelated, we propose the following ansatz : 

4>(x) = (p + + 5(fi(x) where 5<fi(x) = (fi_f(x) + 4> q g(x) cos(q*x) (18) 

Far from the interface, in the disordered region, f(x) = 1 and g(x) = 0, whereas in the 
lamellar phase region, f(x) = —1 and g(x) = 2. These are the boundary conditions for 
the differential equations of the interfacial profile. The slowly spatially varying functions 
f(x) and g(x) play the role of the two order parameters. The difference in the scaling 
of 0_ and <f) q with e close to the tricritical point, enables us to decouple the behavior of 
these two functions. 

Expanding the chemical potential up to the third order in </>_, we find : 

/I - t(j)+ - (j) 3 + = Q!0_ + 30+0 2 + 3 (19) 

where a is define as: 

d 2 F n , 2 3 27 2 e 3 

a = lt = t + 3 ^ = 1 -5 £+ l00 £ "6 (20) 
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Incorporating ansatz fllSl) into G L = F L — p,cj) L , we get : 

a 1 
F L - v4> L = + -(5<pf + + - A m* - + 30+0 2 „ + <PI)5<P + (21) 

2{/>- + W - (<f) 2 <?) cos(g* x) - 2gV sin(g* x)] 2 - 2(/>_ + ^(g' cos(q*x) - qg sm(q*x))} 2 

Close to the tricritical point, the two different length scales behave differently. The 
profile is characterized by the width of the interface, £ , which should diverge at t — t tcp . 
However, in the same limit, the modulation wavelength, 2n/q*, remains unchanged, as 
can be checked from the numerical calculations (see Fig. 5). Therefore, taking the average 
of F L — fi(p L over one period, yields : 

(F L -^ L ) = 

G(<f> + ) + 20 2 _[(f ) 2 - iff] + <P%g" - \g? + g' 2 - \g 2 ] + f ^ f + + (22) 

The variational principle with respect to f(x) leads to the relation : 

af + 30 + 0_/ 2 + l^+g 2 + l^fg 2 + <P 2 -f = <* + 30+0- (23) 
Using expressions fll7D for <f>+, 0_ and <p q , we get up to the first order in e: 

m = i - \g 2 + - ^ ( 24 ) 

This relation enables us to find, when applying the variational principle with respect to 
g(x), the following differential equation : 

V = ^g + 3/#_0 + + \<P 2 J 2 g + l<P 2 q g 3 (25) 

which can be rewritten at the lowest order in e as : 

f </'(*) = (g-g 3 + ^-g 5 )e 2 (26) 
6 lb 

Taking into account the boundary conditions for g(x), the analytical solution of this 
differential equation is found to be : 



where U(x) = \J |(a; — xq) e, and xo can be set to zero without lost of generality. Using 
Eq. (p4|) we can now obtain for / 



f(U) = tanh(C/) - A cosh- 2 ([/) e (28) 

The two functions f(x) and g(x) determine the interfacial profile (j)(x) close to the 
tricritical point and are plotted in Fig. 3 and Fig. 4. The width of the interface scales as 
£ ~ e -1 . Namely, the correlation length exponent v = 1 is recovered, in agreement with 
the scaling laws associated with the tricritical point ||. We note that the main reason it 
is possible to decouple the differential equations for / and g (to lowest order in e) is due 
to the different scaling of <p q and 0_ with respect to e close to the tricritical point [see 
Eq. m 



IV. THE DISORDER-LAMELLAR INTERFACIAL TENSION 

The interfacial tension 7, is defined as the excess of free energy due to the interface. 
It can be written as : 
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J T(x)dx — J[G — Go ^ GL }dx = j[F- ~ '~ ~ ~ '~ - 12 5<J)(x)}dx (29) 



Since up to second order in e, the even part of T(x) vanishes, we have to keep terms 
in T up to third order in e. 



- 1) - ^-(« + 30+0-)/ + f 2 4f - 1) + 0-0 + / 3 (30) 



This expression can be written as function of f(x) only : 

T = (1 - f)( A - |0 2 + - - 3 -^ q ) + odd terms (31) 

and to third order in e 

T = (^cosh- 2 (f/) + odd terms)e 3 + 0(e 4 ) (32) 

where the odd terms in the above equations are odd with respect to x — > —x and will 
vanish upon integration over x, yielding 
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7 




(33) 



This result corresponds to the expected tricritical exponent fi = 2, in agreement with the 
hyperscaling relation \x = (d — l)u, as the upper critical dimension for the tricritical point 
is d = 3 [|. 

We note that in order to have the full expression of f(x) up to first order in e, one 
needs to expand g(x) up to first order in e. Indeed, if g(x) = g (x) + gi(x)e, then f(x) 
can be written as : 



Nevertheless, only the leading order terms of g(x) and f{x) are required to compute the 
interfacial tension. 

Numerical studies of this one-dimensional interface without assuming single-mode 
dominance give results in good agreement with the analytical calculation presented above. 
We have discretized the order parameter <p and chosen the grid spacing to be much smaller 
than 1/q* (typically 180 grid points per period). The free energy (|J) can then be mini- 
mized using the conjugate gradient method. 

The free energies for the disordered (isotropic) and lamellar phases (without using the 
single mode approximation) are computed in order to determine the chemical potential for 
the coexisting phases. Then, we compute the free energy associated with the interface, in 
an interval much bigger than the expected interface length, in order to neglected boundary 
effect. The boundary conditions are taken as reflecting (Dirichlet) boundary conditions. 
The same procedure of Ref. @ using Eq. (p9|) is employed to calculate the interfacial 
energies. For small e = t tcp — t, the optimum wave vector is still q* — 1/ v^2, within a 
good precision. The profiles of the interface for e = t tcp — t =0.35 and e =0.03 are shown 
in Figs. 5 and 6, respectively. 

The interfacial tension is plotted on Fig. 6 as a function of t tcp — t. The log-log plot 
in 6a as well as the 7 1 / 2 plot in 6b manifests clearly that our analytical result are: (i) 
in good agreement with the numerical results obtained from the full minimization of the 
interfacial free energy without assuming the single mode dominance; and (ii) that indeed 
we obtained the expected scaling behavior of the tricritical point. 



In this paper we have derived analytical expressions for the interface profiles between 
a simple modulated phase with a lamellar symmetry and a disordered one. Our results 
are valid close to the tricritical point since we assume that the modulated phase is char- 
acterized by its dominant g-mode in the bulk. This assumption is found to be compatible 
with the scaling laws close to the tricritical point. Our analytical (but approximated) 



f{x) = l--g 2 + —g 2 { y - - l)e - g {x) 9l {x) = f (x) + f x (x)e. 



(34) 



V. DISCUSSION AND CONCLUSIONS 
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results are in accord with numerical minimization of the full free energy functional having 
the same one-dimensional symmetry. 

Unlike simple interfaces between two coexisting homogeneous phases, here the interface 
depends on two separate length scales. Close to the tri critical point, the interface width 
£ diverges and hence gets decoupled from the periodicity of the lamellar phase. It is 
interesting to note that in this case, the interfacial region can involve quite a number of 
lamellae periods as is shown in Figs. 4 and 5. For e =0.03 the interfacial region spans 
over 100 lamellae. 

Our results should be compared with the ones obtained by Fredrickson and Binder 
|| since they obtained a similar envelope equation for the profile (eq. (27) for g in our 
notation). In Ref. the metastable interface close to a symmetric critical point is 
considered. Due to the symmetry, they have / = (in our notation). The high order 
terms (0 5 ) in their profile equation arises from the first-loop correction (order 6 ) in their 
mean-field effective Hamiltonian. In our case, due to the asymmetry, the coupling of g 
to / generates the high order terms even within mean field theory. Needless to say 
that the scaling of the correlation length and surface tension is different in the two models 
because of the difference between critical and tricritical behavior. 

It will be interesting to try to obtain experimental data on lamellar to disorder inter- 
faces close to the tricritical point. Unfortunately, in many experimental realizations the 
one-dimensional phase diagram studied here is modified. The lamellar to disorder tran- 
sition becomes metastable since it is pre-empted close to the tricritical point by phases 
with hexagonal and cubic symmetry 0. Then, the lamellar-to-disordered interfaces can 
exist in thermodynamical equilibrium but only far from the tricritical point. Close to the 
tricritical point the lamellar to disordered transition describes only a metastable situation 
which is of importance, for example, for systems with slow dynamics like copolymer mix- 
tures We hope the systematic expansion we present can be applied to other more 
complicated situations where it is also possible to decouple the important length scales in 
the system. 
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FIGURE CAPTIONS 



Fig. 1 The calculated phase diagram from the modulated phase free energy F L 
as function of the reduced temperature - average order parameter (t, (fi ). The phase 
diagram is symmetric with respect to O — > — <p . The two symmetric tricritical points are 
denoted by tcp. Below the tcp temperature, there are two regions of coexistence of the 
lamellar phase and disordered phase (L+Di) and (L+D2). At even lower temperatures, 
the lamellar phase disappears. We then recover the usual coexistence region between two 
different disordered phases: Di+D 2 ||. 

Fig. 2 Same as in Fig. 1 but as function of the chemical potential - reduced temperature 
t). The continuous line is a second-order transition line. The dashed line is a first- 
order one. They meet at two tricritical points (tcp). The two transition lines separate the 
lamellar phase region from the disordered region. The phase diagram is symmetric with 
respect to /i — > —fi. The /1 = first order transition line meets the two other first order 
lines at a triple point where all three phases: L, Di and D2 coexist. 

Fig. 3 (a) Analytical profiles of the two functions f(x) and g[x) for e = t tcp — t— 0.35. 
The plot scans the entire interfacial region between a bulk lamellar phase (left side where 
g — > 2, / — > — 1) and the disordered phase (right side where g — ► 0, / — ■> 1). (b) Analytical 
interfacial profile <p(x) using the same functions f(x) and g(x) of (a). The value of the 
average concentration in the two coexisting phases is: 0o — 0.25 for the lamellar phase 
and 4>o = 0.50 for the disordered one. 

Fig. 4 Numerical interfacial profile for e = 0.35 and \x = 0.322. In order to demonstrate 
that the most dominant bulk mode remains unchanged even in the interfacial region, the 
bulk lamellar phase is plotted as a dotted line. The horizontal line just below 0o — 0.5 
denotes the average concentration of the bulk disordered phase. 

Fig. 5 Numerical interfacial profile for e = 0.03 and fi = 0.249. The horizontal line 
just above O — 0.3 denotes the average concentration of the bulk disordered phase. The 
width of the interface is much larger than 2n/q*, and over 100 lamellae are located in this 
interfacial region. 

Fig. 6 (a) Interfacial tension as function of e = t tcp — t plotted on a log- log plot. 
The straight line is the analytical interfacial tension 7. The dots represent the numerical 
results as is discussed in the text. The linear dependence of both the analytical and 
numerical results are in accord with the general scaling behavior expected close to the 
tricritical point, (b) Same results as in (a) but plotted differently. The square root of 
the interfacial tension 7 is plotted as a function of e — t tcp — t. The linear dependence is 
another confirmation of the tricritical point scaling 7 ~ e 2 
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